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INTRODUCTION 


Consider  an  inventory  system  using  e continuous  review  (s,S)  policy 
with  constant  lead  time  of  size  T.  The  time  between  successive  demands 
is  ldd  with  distribution  function  H(>)  and  pdf  h(*) • Likewise,  demand 
sizes  are  lid  Integer  valued  random  variables  with  probability  function 
b(«).  All  demands  are  backlogged  until  filled.  We  will  derive  the  dis- 
tribution of  customer  waiting  time.  Since  It  Is  reasonable,  we  take 
s >_ - 1,  which  means  that  no  customer  will  ever  wait  more  than  T. 

While  the  (s,S)  continuous  review  Inventory  system  has  been  greatly 
studied,  there  has  been  little  work  on  customer  waiting  time.  Sherbrooke 
[3]  derived  the  waiting  time  distribution  for  the  special  case  (S-1,S) 
system  subject  to  compound  Poisson  demands.  Simon  [4]  derived  an  expres- 
sion for  the  expected  wait  In  the  (s,S)  system  when  the  demand  process  is 
simple  Poisson  which  can  be  shown  to  be  a particularization  of  L - XW. 

The  most  general  analysis  of  the  (s,S)  continuous  review  system  has 
been  done  by  Sahin  [2]  who  developed  expressions  for  both  the  time  de- 
pendent and  stationary  distributions  of  net  Inventory,  i.e.  on  hand  minus 
backorders,  and  inventory  position,  i.e.  net  Inventory  plus  on  order,  using 
a renewal  - theoretic  structure.  His  advancement  over  earlier  work  was 
In  permitting  demand  size  to  be  a random  variable.  Urbach  [6]  also  analyzed 
a similar  system  under  the  condition  that  no  more  than  one  order  is  out- 
standing. He  was  Interested  In  the  case  of  random  lead  times. 

For  the  sake  of  presentation,  we  first  develop  the  waiting  time  dis- 
tribution for  the  case  of  unit  demand  size.  We  then  use  the  logic  and  some 
of  the  results  derived  to  extend  to  the  case  of  random  demand  size. 


any  assets  ordered  after  y will  not  ta  available  until  after  y + T . This 
means  thaf  a customer  who  arrives  at  t will  wait  <_  * iff  he  receives  one 
of  the  assets  on  account  it  t + i - T,  The  customer  will  get  an  item 
from  A(t+t  -T)  only  if  the  previous  demands  for  those  assets,  d(c+x  -T,t), 
are  less  than  A(t+x  -T) . We  have  then  that 


Pr[W(t)  <x]  - E Pr[W(t)  < x|A(t+x-T)  - a]  Pr[A(t+x-T)  - a] 
a-R+1 


■ £ Pr[d(t+x-T,t)  < a|A(t+  -T)  - a.  Demand  at  t] 
a-R+1 

Pr[A(t+x-T)  - a] 


- £ Pr[A(t+f-T)  - a,  d(t+f-T,t)  < a|Demand  at  t] 
a-R+1 


For  finite  time,  A(t+x-T)  and  d(t+x-T,t)  may  be  dependent  random  variables 
since  knowledge  of  A(t+x-T)  nay  provide  information  about  the  demands  from 
the  start  of  the  inventory  system  until  t+x-T  which,  in  turn,  affect  the 
likelihood  of  d(t+x-T,t).  However,  in  Appendix  A we  show  that  A(t+x-T) 
and  d(t+x-T,t)  are  independent  in  the  steady  state, and  with  the  given 
condition  of  a demand  at  t have  probability  functions 

Pr[A-a]  - ^ ; a - R+l,  R+2....R+Q 


lim  Pr[d(t+x-T ,t)<  d | demand  at  t]  - l-H.(T-x) 

I »oo  <3 


In  other  words,  the  demands  in  the  T-x  units  preceeding  the  present 
customers  arrival  form  an  ordinary  renewal  process  in  the  steady  state. 
Letting  F (•)  denote  the  steady  state  distribution  of  waiting  time  we  then 
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have  from  (1)  that 


«>  V>  ■ 5 ^ - \+k<I-')l  • o i 


T < T 


F <T)  - 1 
w 


As  a matter  of  interest  a special  case  of  the  results  in  Appendix  B 
is  that  E(W)  ■ E(B)/X  where 

E(W)  - expected  waiting  time 
E(B)  * expected  steady  state  backorders 
and  1/X  ■ expected  time  between  demands.  Of  course,  this  is  simply 
an  example  of  L ■ XW. 

Waiting  Time  for  Random  Demand  Size 

For  the  unit  demand  size  case  the  meaning  of  waiting  time  was  obvious. 
In  extending  to  random  demand  size,  we  have  the  problem  of  defining  customer 
wait.  For  example,  what  is  the  wait  when  a customer  who  demanded  10  units 
receives  five  units  immediately,  but  waits,  say  10  days,  before  receiving 
the  other  five  units?  Recognizing  that  the  definition  of  waiting  time 
should  depend  upon  the  context  in  which  the  statistic  is  to  be  used,  we 
avoid  the  problem  of  defining  waiting  time  by  deriving  the  distribution 
of  wait  separately  for  each  unit  in  the  demand.  Later  we  show  the  richness 
of  this  approach  by  demonstrating  how  this  distribution  can  be  used  to 
develop  several  common  performance  measures. 

As  before  we  take  a demand  arrival  to  occur  at  t,  but  allow  the  demand 
size  D to  be  > 1.  Each  unit  in  the  demand  is  identified  by  an  index  j from 
1 to  U.  The  unit  will  wait  <_  t iff  the  demands  preceding  the  unit 
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which  are  vying  for  A(t+v-T)  are  leas  than  A(t+T-T) . In  this  case,  those 
demands  are  the  j-1  units  of  the  present  demand  plus  d(t+t-T ,t) . So 


th  R+Q 

(3)  Pr[j  unit  waits  <_  t]  - E Pr[d(t+r-T,t)  <_  e-J  |A(t+r-T  - a,  Demand  at  t] 

a-R+1 

• Pr[A(t+f-T)  - a] 

R+Q 

■ E Pr[A(t+T-T)  - a,  d(t+T-T,t)  < a-j |Demand  at  tl 
a-R+1 


Again,  as  is  shown  in  Appendix  A,  A(t+-r-T)  and  d(t+r-T,t)  are  independent 
in  the  steady  state.  Several  authora,  [2],  [6],  and  [8],  have  shown  how 
to  compute  the  steady  state  distribution  of  assets.  As  in  the  unit  demand 
size  case,  the  number  of  demand  occurrences  in  the  preceding  T-t  units  form 
an  ordinary  renewal  process  in  the  steady  state.  That  is 


p(d)  - lim  Pr[d(t+x-T,t)  - d|Demand  at  t] 

t-H» 


d 

“ E b (d)[H  (T-t)  - H ..(T-t)];  d > 1 
n-1  n n 0+1 


and 

p(0)  - lim  Pr[d(t+T-T,t)  - OjDemand  at  t]  - 1 - H(T-t) 
t-x» 

Lettlng  jRw('*)  denote  the  steady  state  distribution  function  of  waiting 
time  for  the  Jth  unit  we  have  from  (3)  that 


R+q 

(4)  " 1 Pr[A-a]  P(a-J)  ; 0 < T < T 

3 a-R+1 


J (T)  - 1 

J » 

where  P(a-J)  - E p(k)  if  a - J > 0 

lc-0 

- 0 otherwise 
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Waiting  Time  and  Some  Common  Inventory  Measures 


Expected  number  of  units  backordered  and  initial  fill,  l.e.  the 


fraction  of  demand  satisfied  without  backorder,  are  two  commonly  used 


Inventory  measures.  Waiting  time  relates  to  each  of  these  in  a similar 


way.  We  show  in  Appendix  C that 


Xl-.F  (T))dT 


Expected  Units  Backordered 


Initial  Fill  - I (• 

3-1 


Both  have  the  common  term 


which  has  a simple  Interpretation 


Suppose  a series  of  N demands  is  observed  and  the  units  in  each  of  the 


demands  are  Indexed  as  before.  Let  n. (N)  be  the  number  of  demands  of 


size  i in  the  N demands.  Then  $(N,J)  • E n. (N)/E  i n (N)  is  the  fraction 

i-J  1 i-1 

of  total  units  in  the  N demands  which  have  index  j . In  the  limit  as  N goes 


to  infinity  we  get 


Z b(i)/0 

1-3 


which  is  just  the  proportion  of  units  demanded  which  have  index  j.  Thus 


both  expected  units  backordered  and  initial  fill  are  equivalent  to  taking 


the  corresponding  measures  for  each  possible  unit  in  a demand,  weighting 
by  the  proportion  of  times  that  unit  occurs,  and  averaging. 

There  is  no  particular  advantage  to  actually  computing  the  above 
measures  using  jFw(*).  However,  this  does  indicate  how  simply  some 
measures  are  able  to  be  expressed  with  Fw(«).  Other  examples  are  the 
probability  that  a demand  is  completely  filled  without  waiting  and  the 
expected  number  of  demands  backordered  when  a demand  is  counted  as  back- 
ordered until  completely  filled. 


APPENDIX  A 


DISTRIBUTION  OF  A(t)  AND  d(t,t+z)  GIVEN  A DEMAND  AT  t + l 

The  R,R+Q  Inventory  system  is  observed  from  timeO  to  t+z.  Arbitrarily, 
we  take  the  demand  process  to  begin  at  time  0 with  A(0)  - R+Q.  We  will 
find 

lim  [Pr[A(t)  - a,  d(t,t+z)  - d|A(0)  ■ R+Q,  Demand  at  t+z] 
t-*» 

The  expression  "Demand  at  t+z"  is  used  to  mean  that  a demand  occurs  within 
dz  of  t+z.  Let  N(0,t)  be  the  number  of  demand  occurrences  in  (0,t)  and 
N(t,t+z)  be  the  number  of  demand  occurrences  in  t , t+z . Given  the  values 
of  N(0,t),  N(t,t+z),  and  A(0) , then  A(t)  and  d(t,t+z)  are  uniquely  deter- 
mined by  the  sequence  of  demand  sizes  associated  with  those  demand 
occurrences.  If  N(0,t)  ■ m,  there  is  a countable  though,  in  general, 
infinite  number  of  demand  size  sequences  which  will  result  in  A(t)  ■ a 


given  A(0)  ■ R+Q.  Call  X(m,a,R+Q)  the  set  of  all  such  sequences.  And, 
if  N(t,t+z)  - n then  d(t,t+z)  ■ d iff  the  n demand  sizes  sum  to  d.  We 
nave  then 


(Al)  G(a,d,z,t)  ■ Pr[A(t)  ■ a,  d(t,t+z)  ■ d|A(0)  - R+Q,  Demand  at  t+z] 

d CD 

■ I E Pr[A(t)  “ a,  d(t,t+z)  • d|A(0)  - R+Q,  Demand  at  t+z, 
n*0  ufO 

N(0,t)  » m,  N(t,t+z)  ■ n] 


E b Ad)  £ 
n-0  m»0 


Pr[N(0,t)  - m,  N(t,t+z)  - n|A(0)  - R+Q,  Demand  at  t+z] 

Pr[X(m,a.R+Q)l  Pr[N(o.t)  - m.  N(t,t+z)  - n.  Demand  at  t+z] 

r(t+z)dz 


- F(a,d,z,t,R+Q)  dz/  r(t+z)dz 


where  r(t+z)  - pdf  of  a demand  occurrence  at  t+z 
(also  called  the  renewal  denalty) 

and  Pr[X(m,a,R+Q)]  ■ probability 

a demand  alee  sequence  from  the  eat  X(m,a,R+Q)  occurs. 

Consider  P(m,n,z,t)  - Pr[N(0,t)  - m,  N(t,t+z)  - n.  Demand  at  t+z] 
For  ■ > 1 we  have 

t 

P(m,n,z,t)  - / P(m-l,n,z,t-y)  H(y)dy 
y-o 

and  on  taking  Laplace  transforms  we  get 
P(m,n,z,s)  - P(m-l,n,z,s)h(s) 

Applying  this  recursively  from  m • 1 yields 

(A2)  P(m,n,z,s)  - P(0,n,z,s)[h(s)]“  ; m >_  1 

Consequently,  the  Laplace  transform  of  F(a,d,z,t)  is 

- d - • _ 

F(a,d ,z,s)  - Z bn(d)  P(o,n,z,s)I  Pr[X(m,a,R+Q) ] [h(s) ]“ 
n-o  m-o 

- F^s)  F2(s) 

where  . 

~ a 

Fj^Cs)  - Z bn(d)  P(o,n,z,s) 
n-o 

end  “> 

F2(e)  - Z Pr[X(m,a,R+Q)] [h(s) ]m 
m-o 

Now 

11m  G(a,d,z,t)  ■ 11m  F(a,d,z,t)dz/llm  r(t+z)dz 
t**  f*-> 

•F.(s) 

“ Ha  <— J ) 11m  (sF.  (a) /11m  r(t+z) 

e-*t>  »+o  t-*» 


r,  d - liM  1 Pr[X(m,a,R+Q)h  (t) 

• [11m  / Z b (d)P(o,n,z,y)dy][  t-»-  m-o * ] 


t-*»  o n-o 
provided  the  separate  limits  exist 


11m  r(t+z) 
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Consider  the  first  term 


t d 


11m  / E b (d)P(o,n,z,y)dy  - E b (d)llm  / P(e,n,z,y). 
o n«o  n“o  n t o 


t t z 

But  11m  / P(o,n,z,y)dy  - 11m  / / h (x)h(y+z-x)dxdy 

t-*«  o t-*»  y-o  x»o  n 

z t+z-x 

- / h (x)  11m  / h(y)dy  dx 

x-o  t-*»  y»z-x 


hQ(x) (l-H(z-x))  dx 


- h*  (1-H)  - Hq(z)  - H^z) 


which  Is  just  the  probability  that  the  number  of  renewals  In  an  Interval 
z for  an  ordinary  renewal  process  equals  n. 

Now  consider  the  second  term* 


11m  £ Pr[X(m,a,R+Q)][h  (t)l/lim  r(t+z) 

t-*»  m-o  t-H» 

lia  Pdf  of  a transition  to  asset  state  a at  t 
pdf  of  a demand  occurrence  at  t+z 

t*”° 

-I  /I.Jl 
w.  w y. 

where 

” average  time  between  transitions  to  state  a 
and  y - average  time  between  demands. 

and  which  follows  because  the  number  of  transitions  to  state  a 
is  Itself  a renewal  process.  Moreover,  since  y ■ expected  amount  of  time 


11 


I 


the  system  spends  in  state  a between  transitions  to  state  a,  then 


y/pa  - lim  Pr[A(t)  - a]  - Pr[A  - a] 
t-*» 

Summarizing  then,  we  have  shown  that 


lim  G(a,d,z,t)  - Pr[A-a]  l b (d)  [H  (z)  - H , (z)  ] 

t-H»  n.Q  n n n+1 

This  contrasts  to  Sahln's  [2]  result  for  the  steady  state  distribution 
of  A(t)  and  d(t,t+z)  without  the  condition  of  a demand  at  t+z.  This  con- 
dition changes  the  demand  process  in  a interval  of  length  z from  a 
equilibrium  demand  renewal  process  to  an  ordinary  demand  renewal  process. 


APPENDIX  B 


RELATIONSHIP  OF  WAITING  TIME  TO  EXPECTED  UNITS  BACKORDERED  AND 

INITIAL  FILL 

From  Sailin' 8 work  [2]  we  have  for  the  stationary  system  that 


(BI)  Expected  units  backordered  - E[B  ] * I Pr [Units  Backordered  > j] 

J-o 

as  R+Q  • 

- E E Pr[A-a]  E [H’(T)  - H'  (T)]  C <a+j) 
j-o  a-r+1  n-1  n 0+1  n 


where  H'(T)  -±  / (l-H(y))  H -(T-yJdy 

n u n-i 


and  C (k)  - 1 - E b (i)  “ 1-Un(k)  * n 1 k 

i-1 


- 1;  otherwise 


Letting  Z (T)  - E[B^]  and  taking  Laplace  transforms  we  get  after  some 


algebra  that 


0 


R+Q 


(B2)  Z(S)  - ^ E (l-H(a))  E + E (H(s)n  - H^)11*1) 


a-R+1 


j-o 


n-1 


E tcn+i(**’J)  ~ cn(a+J)I 
0 


where  we  have  used  that  the  Laplace  transform  of  H^(t)  is  (s) 


n-1 


Returning  to  the  time  domain  we  get 
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+ l (H(T-t)  - H(T-t)) 
n-1  n n+1 


[ (I-H(T-t))  e 


f Cn+1<«+J>“  Cn(a+J)) 


It  is  possible  to  show  by  algebraic  manipulation  of  (B3)  that 

It.].;  i=*a=a-  afd-jVx)* 

u J-l  e “ 0 


but  It  is  more  appealing  to  argue  the  above  result  straight  from  the  terms 

c.+i(*+J)  " C (e+j) 

In  (B3).  The  term  Z Is  the  proportion  of  units  in  the 

st  J“° 

n+1  demand  which  wait  longer  than  t given  that  A(t+r-T)  - a,  and  that 


the  number  of  demand  occurrences  in  (t+r-T ,t)  - n.  In  this  perspective 


the  n+1  demand  Is  the  customer  arriving  at  t.  By  probablistlcally 


weighting  over  all  values  of  n,  we  get  the  proportion  of  units  in  the 


arriving  customers  demand  which  must  wait  longer  than  t.  But 
« T 

E l-B(l-l)  / (l-.F  (t) )dx » by  the  arguments  given  in  the  report  Is  the 
j-l  9 0 J w 

same  proportion. 


Initial  fill,  IF,  is  tha  fraction  of  total  demand  which  Is  filled 


without  wait . Thus 


Number  of  units  filled  without  wait 
Number  of  units  demanded 


N(t)  - number  of  demands  In  period  t 


amount  of  1 demand  filled  without  welt;  1 - 1,  2 


So 


N(t)  N(t) 

IF  - 11m  Z 1 / 11m  Z d./N(t)  . 

t-H»  1«1  N(t)  t-H»  1-1  1 


. Expected  number  of  units  filled  without  wait  per 

Q ■ ' with 

probability  1. 

Now  the  expected  number  of  unite  filled  without  wait  per  demand  Is 

j«l  tPr?  that  J °r  “°re  Unlt*  ln  * deoand  are  Filled  without  wait]  which 
equals  (jFn<0)) (l-B(j-l)) . In  other  words,  for  j or  more  units  to  be 
filled  Immediately,  there  must  be  at  least  j unite  demanded  and  at  least 
the  J unit  does  not  wait.  Consequently 

if  - e U-BU-DX^CO)) 

3-1  0 
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